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The classical and quantum dynamics for an ra-dimensional generalization of the kicked planar 
(n = 1) rotator is considered. The classical behavior of this system is, in essence, described by a 
one-dimensional kicked rotator in an additional effective centrifugal potential. Therefore, typical 
phenomena like the diffusion in classical phase space are similar to that of the one-dimensional 
model. For the quantum dynamics such a result is not expected as in this case the evolution does 
depend in a very complicated way on the number n of degrees of freedom. In the limit n — > oo we 
find the free undisturbed quantum motion. For finite values of n (I < n < 26) we study numerically 
the quantum dynamics. Here, we always find localization independent of the actual number of 
degrees of freedom. 
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I. INTRODUCTION 

m 

ON 

Quantum signatures of classically chaotic systems have 
i-Hattracted increasing attention in the last 10-15 years 
^[ 1-5] . One particular class of systems which has been 
^investigated is characterized by time-periodic Hamilto- 
^Qnians [6]. The most prominent member of this class is 
jthe so-called kicked planar rotator. This is a point parti- 
|de which moves freely along the unit circle in the plane. 
This motion is perturbed by time-periodic kicks whose 
O^strength does depend on the actual position of the parti- 
^cle. Accordingly, the Lagrangian of this system may be 
(^written as 
OO 

o 
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fc^Here the angle $ £ [0,27r) describes the position of the 
i-^particle and I stands for its momentum of inertia with 
I respect to the origin. Furthermore, r > is the pe- 
Oriod of the kicks and K > is the maximal kicking 
^ strength. It is this "stochasticity" parameter K which 
^controls the classical dynamical behavior of the system. 
Near Kt 2 / 1 k, 0.9716 . . . there is a transition from local 
to global stochasticity taking place [7]. This is typically 
visible in the linear increase of the system's energy in 
time. 

In the present paper we will introduce an n- 
dimensional generalization of the system characterized 
by (1) and study its classical as well as its quantum dy- 
namics. More precisely, we consider the free motion of a 
point mass on the n-dimensional unit sphere embedded in 
the (n + l)-dimensional Euclidean space. This motion is 



periodically disturbed by kicks, whose strength depends 
only on one of the n + 1 Cartesian coordinates. On the 
classical level this generalized system is essentially identi- 
cal to the one-dimensional the effective equation 
of motion does not depend on the number of degrees of 
freedom. This is in contrast to the quantum dynamics 
of the n-dimensional kicked rotator, which explicitly de- 
pends in a rather complicated way on the number n. The 
study of this n-dependence has been our primary motiva- 
tion for this work. Let us also mention that our model for 
n = 2 is somehow related to the kicked top extensively 
discussed by Haake and coworkers [8]. See also ref. [5]. 

Our paper is organized as follows. In the next Section 
we explicitly define the n-dimensional generalization of 
(1) which we are going to investigate. Section 3 is de- 
voted to the study of the classical dynamics. We inte- 
grate the classical equation of motion and obtain a cer- 
tain n-dependent generalization of the well-known stan- 
dard map. The observed diffusive behavior of this map 
is compared with that for the one-dimensional case. In 
Section 4 we will present analytical as well as numerical 
results for the quantum dynamics. In the limit n — ► oo 
the n-dimensional kicked rotator behaves like a free ro- 
tator (without kicks). For finite values of n we always 
observe — by numerical methods — the phenomenon 
of so-called localization. The corresponding localization 
length is, in contrast to what one might expect, found to 
be independent of n up to values as large as 26. Unfor- 
tunately, our numerical method does not allow to draw 
firm conclusions for even larger values of n. 
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II. THE MODEL 

The n-dimensional generalization of the planar rota- 
tor is simply the free motion of a point mass on the unit 
sphere embedded in the (n + l)-dimensional Euclidean 
space 



(Xl 



) gir 1 



n + l 



n + 1 



(2) 



The unit sphere S n is usually parametrized by polar co- 
ordinates: 



x r = e r sin r 

Z = X n + l = cos 



1,2, . . .,n, 

< < 7T. 



(3) 



Here $ denotes the angle between the unit vector x := 
(x\, . . . , x n+ i) and the z-axis which is chosen to be the 
coordinate axis for the last component of x. In the above 
e; stands for the i-th component of an arbitrary unit 
vector e := (ei, . . . , e n , 0) G IR n+1 perpendicular to the 

z-axis. Obviously, x = i? 2 + e sin 2 $ is the squared 
modulus of the velocity of a point particle moving along 
a path x = x(t) G S n . 

The Lagrangian of the n-dimensional kicked rotator is 
then defined by 



L n (d,d,e,t) 



K cos d 



■_'2 

e sin 



2 <? 



(4) 



where M := \M\. Let us note that for the special case 
M = the above Lagrangian is identical to that for the 
planar rotator. Indeed, classically the case M = can- 
not be distinguished from the system (1) as in the former 
case the motion takes place in a plane containing the z- 
axis. Let us add, that in this case the parametrization 
(3) is singular and should be replaced by < d < 2tt. 
In the general case M > the Lagrangian (6) may be 
interpreted as that for the kicked planar rotator in an 
additional time-independent (effective) potential. 

The classical equation of motion, which may be de- 
rived form L e f f , can easily be reduced to a discrete-time 
map on the two-dimensional classical phase space SI := 
[0, 7r] x HI. This is because (6) describes a regular motion 
between two successive kicks. Introducing the quanti- 
ties Zj := z(jr) = cos^(jr), the position (z-coordinate) 
at the j-th kick, and Zj := lim £ |o z(j T + £ )> the- corre- 
sponding velocity just after the j-th kick, the equation 
of motion leads to the phase-space map ($j := arccos Zj, 
dj := -Zj/smdj) 



T 



(7) 



This map is explicitly given by 

Zj+i = Aj sin [ujjT sgn (ij) + <f>j]. 



LOjAj sgn (ij) cos [ujjT sgn (ij) + <f>j] 



(8) 



-k(1 -zj+i), 
where we have set 



As in (1), I is the moment of inertia, K is the kicking 
strength and r the period of these kicks. Due to the ro- 
tational invariance of (4) about the z-axis we have n — 1 
constants of motion: 



M := dhnjde = Iesm 2 tf. 



(5) 



They allow the reduction of system (4) to an effective 
one-dimensional system, which in turn allows for a so- 
lution of the classical equations of motion in terms of 
a time-periodic classical map in a reduced phase space. 
We will derive this n-dimensional generalization of the 
standard map in the next Section. 



III. THE CLASSICAL DYNAMICS 

Due to the constants of motion (5) the dynamics of the 
system (4) can be reduced to a one-dimensional system 
characterized by the effective Lagrangian 



2 27 sin 2 



-K cos d 



6 I - 



(6) 



M' 2 



/ (M/J) 2 + (z J ) 2 
1-z 2 



arcsin ( ) , sgn(z) 



z/\z\ for z ^ 
1 for z = 



(9) 



and introduced the rescaled kicking strength k := Kt/I. 
Note that 



lim( z(]T + e) 

£|0 V 



z(jr -£))=« (z 2 - 1) 



(10) 



We leave it to the reader to show that for M = the map 
(7), respectively, (8) reduces to the well-known standard 
map for the planar rotator. 

We have iterated the map (7) numerically for various 
values of the parameters k and M using units such that 
I = r = 1. Figures 1-3 show typical phase-space por- 
traits generated by the above map. It is obvious that with 
an increasing value of k, keeping M fixed, the islands of 
regular motion become smaller. Whereas, an increas- 
ing M for fixed k will cause growing islands of stability. 
In particular, Figure 1 displays, because of M = 0, the 



2 



10 



•10 













;;<::::: 3 C. 


























^ 






:'!>c.:;-;3: C" 








7r/2 ^ 7T 



FIG. 1. Phase-space portrait generated by the map (7) 
for parameter k = 1.5 and M = 0. 
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FIG. 2. Same as Figure 1 for k = 1.5 and M = 1. 



the energy £j := JwJ/2. For sufficiently large k it is ex- 
pected that a random-phase approximation is valid [7]. 
That is, one assumes a uniform distribution for dj in 
the interval [0, 7r] and an independent distribution for dj 
which is symmetric about the origin. These assumptions 
directly lead to a linear (hence diffusive) growth of the 
energy with time. In other words, the diffusion coefficient 

D(k) := lim(E j -E )/j (11) 
exists. Its leading dependence on k for large k is given 

by 

D(k) ~ Ik 2 /4 =: D (k), k^oo. (12) 

For the calculation of subleading corrections one has to go 
beyond the random-phase approximation. For the spe- 
cial case M = this has been done by Rechester and 
White [9] (cf. also [7]) with the result 

D(k)/D (k) ~ 1 - 2 (j 2 (/s) + Jf (k) - Jl(/s) - Jl(/s)) , 

(13) 

where J„ denotes the Bessel function with index iA 

We have calculated the k dependence of the diffusion 
coefficient (11) for various values of M with the result 
that D(k) is well described by (13) independent of the 
parameter M . Figures 4 and 5 present our results for 
M = and M = 10, respectively. The observed enhance- 
ment in the diffusion coefficient near k ps 27rm, m G IN, 
is attributed to so-called acceleration modes [10]. 

In concluding this Section we may say that the classi- 
cal n-dimensional kicked rotator is, as far as its stochas- 
tic nature is concerned, identical to the one-dimensional 
kicked rotator. 



tt/2 



FIG. 3. Same as Figure 1 for k = 2.5 and M = 1. 

well-known phase-space portrait of the kicked planar ro- 
tator. 

In order to study the stochastic nature of the map (7) 
quantitatively we have looked at the diffusive growth of 



IV. THE QUANTUM DYNAMICS 

Let us now consider the quantum version of the model 
introduced in Section II. It is characterized by the Hamil- 
tonian: 



21 



A + K cos 6 ^2 s 



(14) 



Here A stands for the Laplace-Beltrami operator defined 
on the unit sphere S n and represents the position op- 
erator corresponding to the polar angle 

Before we study the unitary time evolution generated 
by (14) we comment on the operator A and the Hilbert 
space L 2 (S") it is acting on. It is well known [11] that 
L 2 (S") can uniquely be decomposed into invariant or- 
thogonal subspaces V 1 each of which carries a unitary ir- 
reducible representation of the group SO(n + 1) labeled, 
in analogy to 50(3), by I = 0, 1, 2, . . .: 
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FIG. 4. The diffusion coefficient D(k)/D (k) for M = 0. 
The solid line shows the expected analytical behavior (13). 



i 2 (5") = 0D' 



£=0 



The £-th subspace has the finite dimension 

T(£ + n- 1) 



{21- 



1) 



r(£ + i)i» ; 



(15) 



(16) 



where T denotes Euler's gamma function. In fact, the 
subspace V 1 is the eigenspace of A corresponding to the 
eigenvalue — £(£ + n — 1), see ref. [11]. An orthonor- 
mal basis in V 1 is the common eigenbasis 
of the n — 1 commuting Casimir operators of the sub- 
groups SO(p), p = 2, . . . , n, appearing in the group chain 
50(2) C • • • C SO(n) C SO(n + 1) [12]. Here M stands 
for the (n — l)-tuple 



M := (mi,m 2 , 



m n _ 2 , m n _i 



where 



\mi \ < m 2 < • • • < m n _ 2 < m n _ 



< 



with mi G TL and m r G INo for r = 2, 
spectral decomposition of A now reads 



A 



E 

1=0 



£(£■ 



l)Y\£M)(£M\. 



M 



(17) 

(18) 
The 

(19) 



With the above remarks we can discuss the time evo- 
lution generated by the Hamiltonian (14). Let us first 
consider the one-kick evolution operator 



U 



exp 



iTir 
~2I 



■A 



exp 



Kt cos 



(20) 



which describes the time evolution for one period r start- 
ing right before a kick. If denotes a state right before 
the j-th kick the time evolution of this state is given by 



FIG. 5. Same as Figure 4 for parameter M = 10. 

I*i+i> = UlVj). (21) 

The properties of this quantum map are our main con- 
cern. 

In order to make (21) more explicit we expand the state 
\^j), for each j, into the basis {|£.M}}: 



(22) 



'M 



where 



(23) 



Rewriting (21) in terms of these expansion coefficients 
yields the recurrence relation 



ai M {j + 1) = exp { + n - 1) 



J2 ai< M < (j) {£M\ exp {-(i/Ti)KT cos 0} \£'M') 



(24) 



What remains to be done is the explicit calculation of 
the matrix element appearing in (24). In doing so we 
will work in the coordinate representation 



Y tM (d,e) := (d,e\£M), 



(25) 



that is, with the n-dimensional generalization of the 
spherical harmonics. They are orthonormalized accord- 
ing to 



4 



d$ sin"' 1 tide Y;, M ,($,e)Y iM ($,e) = b tJJ b, 



M,M' 



(26) 



Here and below we will adopt the notation and results 
of ref . [13]. In terms of the spherical harmonics (25) the 
matrix element reads 

(£M\exp {-(i/h)I<T cos 0} \£'M') 

d$ sin" -1 d de exp {-(i/h)Kr cos «?} ( 27 ) 
xY* M ^,e)Y t , M ,^,e). 

In order to make this integral somewhat more explicit 
we expand the exponential function in terms of spherical 
harmonics: 

oo 

exp {-(i/Ti)KT cos = ^ Yt»o(ti, e ) 5^»(A), (28) 

l"=0 

where we have introduced A := Kt/Ti for brevity. Be- 
cause of its invariance under rotations about the z- 
axis only those spherical harmonics with M" = O := 
(0, . . . , 0) contribute to the expansion of this exponen- 
tial. The expansion coefficient is explicitly given by 



Bi(\) 



dd sin"" 1 



$ deYf* ($, e) exp{— iX cos «?} 



n+l 



(_ j y/2 (2/A) ("-l)/2 



x J<+(n-l)/2W- 



(29) 



For the integration we have made use of the fact that 
Yio($,e) is expressible in terms of a Gegenbauer poly- 
nomial. See, for example, ref. [11]. 

In a second step we make use of the equality 



V ' ' 1 IM 

x (£'0; £"0\(£'£")£0)((£'£")£M \£'M';£"0) 



(30) 



which can be derived from eq. (34) in connection with eq. 
(15) of ref. [13]. Here {£ 1 Mi;l 2 M 2 \(h^M) denotes 
a generalized Clebsch-Gordan coefficient and \S n \ := 
27r («+i)/2/ r (n±i) stands for the volume of the unit 
sphere S n . Making use of the orthogonality relation (26) 
we finally arrive at the result 



(£M\exp {-(i/h)I<T cos 0} \£'M') 



xB lll (\)(£'0;£"0\(£'£")£0)((£'£")£M\£'M';£"0). 



(31) 

Using the fact that the Clebsch-Gordan coefficient 
{£'M';£"0\(£'£")£M} vanishes unless M = M' we can 
put the map (24) into the simple form 




ai M {j + 1) = exp { + n-l) 



OO 



J2 <*i'MU)Bi»(\)G M (t',t";t), 



(32) 



l',l" = 

where we have introduced 



G M (£',£";£) 




r —(£'0;£"0\(£'£")£0) ^ 
x((£'£")£M\£'M;£"0). 



This quantity is of purely geometric origin and, in 
essence, controls which angular-momentum eigenstates 
\£M) can be reached from an initial eigenstate |i".M} 
due to an angular-momentum transfer £" stemming from 
the kick. In contrast to this, the term Biji(X) is of dy- 
namical origin and represents the weight with which the 
kick contributes to the angular-momentum transfer. 

The map (32) also explicates that the (n — l)-tuple M 
is conserved. That means, there exists n — 1 constants of 
motion. These are in fact the counterparts of the n — 1 
classical constants of motion (5). In particular, selecting 
an initial state | ) such that (^^(O) oc 8m, o the map 
(32) represents the quantum version of the classical map 
(8) with M = 0, hence, the well-known standard map. 
Whereas (8) does not depend on the number n of degrees 
of freedom, the quantum map (32) explicitly depends in 
a rather complicated way on n. With (32) we have IN dif- 
ferent quantum versions for the classical standard map. 

Let us note that in principal we are now in a position 
where we could study the map (32) numerically. How- 
ever, the explicit calculation for the Clebsch-Gordan co- 
efficients is rather complicated and also requires an exten- 
sive numerical effort [13]. Only the special initial condi- 
tion M = O allows for a practicable numerical iteration 
as in this case the Clebsch-Gordan coefficients are avail- 
able in closed form [13]. Hence, in the following we will 
concentrate on this special case. 



A. The special case M. = O and the large n limit 

In this Section we consider the quantum map (32) with 
the special initial condition M = O . For convenience we 
use the notation ai(j) := aioij)- As we have already 
mentioned, the Clebsch-Gordan coefficients for this case 
are known in closed form [13]. Hence, we are also able to 
express Go(fi,4;4) in a closed form: 



r(j- 



i) 



3 

><n 



r 2 (a±i) r(n - i)r (J + ^f- 
(A- + 2 f i ) T(j-e i + 2=±y 



T( J - £i + 1) 



(34) 
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This expression is only valid if J := (£\ + £2 + £s)/2 6 INo 
and the £i obey the triangular relation £j + £^ > ^ > 
\£j ~ &k\, G {1,2,3}. This conditions are similar 

to those known for the case n = 2. Whenever J is not 
an integer or the triangular relation is not fulfilled the 
Clebsch-Gordan coefficients and hence Go(£-i, £-2', ^3) van- 
ishes. 

With this result the map (32) reduces to 



(35) 



ai(j + l)= exp^-—£(£ + n-l) 



J2 a f ,(j) B f ,,(\)Go(£', £";£)■ 



l',l" = 

This recurrence relation together with the explicit expres- 
sion (34) is now in a form which allows for a numerical 
iteration. 

However, before we are going to report our numeri- 
cal results let us briefly consider the limit n — ► 00. Using 
the asymptotic form of the Bessel function for large index 
and the Stirling formula applicable for large arguments of 
the gamma function one finds the following leading-order 
terms for n — ► 00: 



Biji(X) « y^\d^ (-i\/nf 



r (£' £"-a ~ nl " /2 aA^ + W+i) (36) 

Go{tJ ' J) ~ vm^- iy-*+u (36) 



T( J - £' + 1)T( J - £" + 1) 
Inserting these approximations into (35) results in 
iTir 



a t (j + 1) « exp \ —^j£(£ + n-l) 



\l"=Q 



(-if l 2 (\/y/n) 1 ^/T(£i + l)T(£+ 1) 

r(j — £ + i)r(j -£' + i)r(j - 1" + i) ' 



(37) 



It is obvious that in the limit n — ► 00 the leading term 
corresponds to £" = 0. Due to the triangular relation 
the condition £" = also implies £' = £. Therefore, in 
essence, we end up to leading order in l/n with the free 
evolution: 

a t (j + 1) « exp j-^ £(£ + n - 1) J a,(j). (38) 

This fact may be explained as follows. Whereas in 
the classical limit the motion takes places in a two- 
dimensional plane, say the (x n , « n _|_i)-plane, quantum 
mechanics allows also for fluctuations in the other degrees 
of freedom (x\, . . . , x n -\). Hence, in the large n limit the 
quantum model may be viewed as a planar kicked quan- 
tum rotator coupled to (n — 1) harmonic oscillators. The 
increase, respectively, the decrease in the energy is com- 
pensated by these oscillators. 



B. Numerical results for A4 = O and finite n 

We have numerically iterated the map (35) for vari- 
ous parameter sets and various initial states |\?o)- From 
the resulting angular-momentum distribution |d£(j)| 2 we 
have calculated the expectation value of the kinetic en- 
ergy, 



(^■1-^72^1^ 

t 2 A 



_^( £ + n _l) Mi )|2 



(39) 



1=0 



and the coarse-grained angular-momentum distribution 



Ai(j max ,Aj) 



1 

Aj 



(40) 



In the above A is a cutoff parameter and j max and Aj 
are chosen such that (40) results in a quasi-stationary 
coarse-grained angular-momentum distribution, that is, 
it becomes independent oi j max . As the long-time behav- 
ior of the above quantities is, in essence, independent of 
the initial state we present only results for |M/ ) = |0O). 

As in the one-dimensional case we observe resonances, 
that is 



E(j) <x i 2 



(41) 



if the quantity (Tir/27) is a rational multiple of 27r. A 
typical resonance is shown in Figure 6 where the units 
have been chosen to I = Ti = 1 and r = Ait. Here, 
the angular-momentum distribution |d£(j)| 2 does not be- 
come quasi-stationary. With increasing time higher and 
higher angular-momentum eigenstates are excited. Af- 
ter, say, 500 kicks a cutoff parameter A = 7000 is no 
longer sufficient. The numerical iteration is no longer 
practicable. 

For the non-resonance case we used units Ti = r = I = 
1, whence K = A = k. A small kicking strength, that is, 
K < 4.5, has been found to lead always to quasi-periodic 
oscillations in the energy (39). See Figure 7 for a typi- 
cal example. Only the lowest angular-momentum eigen- 
states are excited and, hence, an angular-momentum cut- 
off A pa 20 is sufficient. 

For a larger kicking strength the energy (39) initially 
increases linearly in time, displaying the classical behav- 
ior (11). After about 20-50 kicks the quantum expecta- 
tion value deviates from the classical one and shows fluc- 
tuations about a constant value (see Figure 8). This phe- 
nomenon is called dynamical localization and is known 
from the one-dimensional quantum rotator. Indeed, the 
corresponding state is localized in (angular) momen- 
tum space. Its coarse-grained angular-momentum dis- 
tribution (40) becomes quasi-stationary, for sufficiently 
large j max and shows an exponential decay in £: 



At ~ exp{a£}. 



(42) 
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FIG. 6. Expectation value (39) of the kinetic energy for a 
resonance. Parameters are ft = 7 = A = 1, r = Air and n = 2. 




FIG. 7. Same as Figure 6 for the quasi-periodic case 
ft=r=7=A=l and n = 2. 




FIG. 8. Same as Figure 6 for parameters ft = r = I = 1, 
A = 9 and n = 2. Here A is sufficiently large to lead to lo- 
calization. The dotted line indicates the classical diffusion 
according to (If). 
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FIG. 9. The coarse-grained angular momentum distribu- 
tion At for parameters ft = r = 7 = 1, A = 14 and n = 2. 
The localization is indicated by the solid line which is a fit to 
the exponential decay (42). 




FIG. 10. The negative inverse localization length a as a 
function of A = k (units ft = r = I = 1) obtained from a fit 
of (42) to our numerical data for n = 2 (0), n = 4 (+) and 
n = 16 (□). The solid line displays the analytical behavior 
(43). 

The constant a is the negative inverse of the so-called 
localization length. A typical coarse-grained angular- 
momentum distribution is shown in Figure 9 for which we 
have chosen j max = 500 and Aj = 50 in (40). The solid 
line is a fit to the exponential behavior (42). According 
to Shepelyansky [14] one expects the localization length 
to be identical with the classical diffusion constant (11), 
that is, 



a = -l/D(n). 



(43) 



Despite the fact that Shepelyansky has derived (43) for 
the n = 1 case we find that this relation also holds for 
n > 1. In Figure 10 we compare the analytical behavior 
(43), where D(k) is given by (13), with our numerical 
results obtained from a fit to (42) for various values of n. 
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The n-independence of the localization length is explic- 
itly visible in Table I where we show numerical results 
for a = a(n) up to n = 26. Within this range of n we 
do not see any signature of a transition to the free mo- 
tion expected from the limit n — ► oo as discussed in the 
previous Section. 

V. CONCLUDING REMARKS 

In this paper we have studied the dynamics of an n- 
dimensional generalization of the kicked planar rotator. 
Classically, this system is found to have the same stochas- 
tic properties, that is diffusion in phase space, as those 
known from the well-studied case n = 1. This is not sur- 
prising as the classical dynamics does not depend on the 
number of degrees of freedom. 

This is in contrast to the quantum dynamics which ex- 
plicitly depends on this number n. Indeed, taking the 
naive limit n — ► oo one recovers to leading order the free 
quantum motion [15]. That is, the kicks are lost in the 
sea of the infinitely many degrees of freedom. However, 
signatures of this disappearance have not been found by 
us for any finite value of n which we have investigated. It 
turned out that we are limited by n < 26 for an explicit 
numerical calculation. In this range for the number of 
degrees of freedom we found no significant dependence 
on n. What we did find are resonances, quasi-periodicity 
and localization known from the quantum dynamics of 
the kicked planar rotator. It is still unclear, whether 
the transition from dynamical localization to free motion 
takes place at a finite value of n or only in the limit 
n — ► oo [16]. So more work is clearly needed for a fuller 
understanding. 
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TABLE I. The negative inverse localization length as a function of the number of degrees of freedom n. Parameters 
are I = h = r = 1 and k = A = 11. Although these numerical values are systematically above the theoretical value 
a = —4.1248 x 10 -2 expected from eq. (43), they do not indicate any ra-dependence of a. 



n 1 


2 


4 


6 


8 10 


12 


14 


16 


18 


20 


22 


26 


-a x 10 2 2.5658 


3.6280 


3.0773 


2.1610 


2.8914 2.3077 


3.1164 


2.9679 


3.2618 


3.3897 


3.0227 


2.5675 


2.9807 
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